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The chiral expansion of nucleon properties such as mass, magnetic moment, and magnetic polarizability 
are investigated in the framework of chiral perturbation theory, with and without the heavy-baryon expansion. 
The analysis makes use of a pion-mass dispersion relation, which is shown to hold in both frameworks. The 
dispersion relation allows an ultraviolet cutoff to be implemented without compromising the symmetries. After 
renormalization, the leading-order heavy-baryon loops demonstrate a stronger dependence on the cutoff scale, 
which results in weakened convergence of the expansion. This conclusion is tested against the recent results of 
lattice QCD simulations for nucleon mass and isovector magnetic moment. In the case of the polarizability, the 
situation is even more dramatic as the heavy-baryon expansion is unable to reproduce large soft contributions to 
this quantity. Clearly, the heavy-baryon expansion is not suitable for every quantity. 
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I. INTRODUCTION 

The dependence of hadron properties on the quark masses 
— the chiral behavior — is crucial for interpreting the modern 
lattice QCD calculations, which usually require an extrapola- 
tion in the quark mass. It is also important for determining 
the quark mass values, as well as for any quantitative descrip- 
tion of chiral symmetry breaking. Chiral perturbation the- 
ory (xPT) (Hi], a low-energy effective field theory (EFT) 
of QCD, should in principle describe the nonanalytic depen- 
dencies on the light-quark masses in a systematic fashion. The 
analytic (series-like) dependencies are more arbitrary in x?T, 
and are specified in terms of low-energy constants (LECs). 
They are to be fixed by matching to the underlying theory, 
which is usually achieved by fitting to lattice QCD results. 

In the baryon sector of xPT one often invokes an additional 
expansion in the inverse baryon masses called the heavy- 
baryon expansion, or HBxPT yllj]- The baryon x?T with- 



out the heavy-baryon expansion will be referred to as BxPT. 
While both expansions should converge to the same result, the 
question is whether they converge in a natural way (cf. |]5l Q 
for an explanation of 'natural'). If the difference between 
them at any finite order is as crucial as is claimed, it is clear 
that they cannot both converge naturally. Therefore, the goal 
is to establish which of HB- or BvPT has the more natural 
expansion. 

Up to a given chiral order, the B%PT result can be written 
as the HBxPT result and a series of contributions, which are 
nominally of higher order in HBxPT. Whether these contribu- 
tions are indeed of smaller size can only be checked in explicit 
calculations, and many cases it has been observed that the dif- 
ference between a given-order of HB- and BxPT results is un- 
naturally large I7l4l4ll. A notable example is provided by the 
magnetic polarizability of the proton, for which the leading 
chiral-loop contribution predicts —1.8 in BxPT B15I1 or +1.3 
in HBxPT [4], in units of 10~ 4 fm 3 . These large differences 
are reconciled in practice by adopting unnaturally large values 
for some of the LECs, appearing at higher orders in HBxPT 
(see e.g. IU6I0 . This is not a solution of course, but rather a 
restatement of the problem. 

In this work, the problem is investigated by using the re- 
cently established pion-mass dispersion relation [ T7J, ll8|] : 
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where the static quantity / is a complex function of the pion 
mass squared rrv\. Throughout this work, the focus is on the 
following static properties of the nucleon: the mass / = Mm, 
anomalous magnetic moment (AMM) / = k^t, and magnetic 
polarizability / = /3jv. The dispersion relation results from 
the observation that chiral loops in x?T are analytic functions 
in the entire complex plane of m\ except for the negative real- 
axis, which contains a branch cut associated with pion produc- 
tion, see Fig.Q] 



for the empirical values of qa — 1-27 and / w ~ 92.4 MeV. 
The nonanalytic term follows from the self-energy graph in 
Fig. 12 which yields the familiar form in HBxPT: 
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FIG. 1. (color online). The complex t — m^ plane, with the branch- 
cut along the negative real axis, and the contour indicating the ana- 
lyticity domain. 



As will be demonstrated in Sec. [II] there is a correspon- 
dence between integrating over t in the above dispersion rela- 
tion and the integration over the momentum in a chiral loop. 
One can cut off the ^-integration without danger of compro- 
mising the symmetries of the theory. Thus an ultraviolet cutoff 
A is introduced in Sec. [HI] After subtracting positive powers 
of A into the available LECs, the aim is to determine the scales 
at which differences between the HB- and B^PT appear. It 
is observed that cutting off the dispersion relation is equiva- 
lent to the 'sharp cutoff version of finite-range regularization 
(FPvR), introduced originally to improve the HB%PT expan- 
sion by resumming the chiral series through the introduction 
of a regulator into the loop integrals, see IU9l - l28ll . The equiv- 
alence between the cutoff dispersion relation and FRR readily 
allows for an extension of the FRR to the realm of B^PT. 

In Sec. lIVI the HB- and B%PT FRR formulae for the cases 
of the nucleon mass, AMMs, and magnetic polarizability are 
obtained at order p 3 , and their residual cutoff-dependence is 
studied. In Sec.[VJ these formulae are confronted with experi- 
mental and lattice QCD results. The recent lattice QCD simu- 
lations of PACS-CS [290 and JLQCD ^ Collaborations are 
used for the nucleon mass, and simulations from QCDSF [31] 
Collaboration are used for the AMMs. The conclusions are 
presented in Sec. [VI] In the Appendix, the dispersion relation 
in the quark mass is considered briefly, and a condition for its 
'subtractions' is discussed. 
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where k is the magnitude of the 3 -momentum running in the 
loop. After a change of integration variable to t = — k 2 , one 
obtains: 
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which is simply the pion-mass dispersion relation of Eq. ([T), 
with: ImS™(i) = \n (— t) 3 ^ 2 . Thus it is not only evident 
that the heavy -baryon loop obeys the pion-mass dispersion re- 
lation, but also that there is a correspondence between the in- 
tegration over the pion mass and the loop momentum. It can 
easily be seen that cutting off the loop momentum at some 
scale A is equivalent to placing the lower-limit of integration 
overt at —A 2 . 

Having observed this equivalence, the dispersion relation 
will be used in preference to the loop integrals themselves. 
The reason for this is twofold. Firstly, in order to express 
the loop integral with only a single /c-integration, the fco and 
angular integrations must be evaluated, which is usually more 
difficult than finding the imaginary part of the loop (especially 
for multi-loop graphs and BxPT expressions). Secondly, it is 
explicitly clear that no symmetries of the theory are compro- 
mised in the cutoff of t. Of course, it has been shown that FRR 
schemes with a sharp cutoff in the momentum are consistent 
with chiral symmetry [25]. However, in ^-integration, consis- 
tency with the symmetries is more explicit, and violations of 
analyticity become apparent immediately. 

The cutoff-dependence of a given quantity is meaningful 
only after the loop contributions have been renormalized. 
When using the dispersion relation, the usual procedure of 
dimensional regularization and cancellation of infinities by 
counter-terms is replaced by 'subtractions'. For instance, the 
nucleon mass requires at least two subtractions at m 2 = 
(cf. Appendix lAl. so that the third-order contribution can be 
written as: 



II. INTEGRATION OVER LOOP MOMENTUM VS. PION 

MASS 

Consider the example of the nucleon mass, which, takes the 
following simple form to order p 3 in HB%PT: 
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Mn and C\ are the LECs to this order, and the coefficient of 
the nonanalytic term is given by: 
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FIG. 2. The pion loop contribution to the self-energy of the nucleon, 
providing the leading nonanalytic contribution to the nucleon mass. 
All charge conserving transitions are implicit. 



with the full result, to order p 3 , given by Eq. (0. The LECs 
play the role of the subtraction constants. 

Introduction of an ultraviolet cutoff in the finite integral af- 
ter the subtractions allows one to separate the low- and high- 
momentum contributions. It can also provide information 
about the scale at which the HB- and B%PT results start to 
deviate. 



of the cutoff any deviation occurs between the HB- and B%PT 
results. If the deviation begins at A -c 1 GeV, then the differ- 
ences between the two expansions cannot be reconciled in a 
natural way. In the next section, this situation is examined us- 
ing several specific examples, and for each of them a different 
picture is obtained (cf. Fig. [3j. 



III. INTRODUCING A CUTOFF: FINITE-RANGE 
REGULARIZATION 

The chiral expansion of an observable quantity / is an 
expansion in the quark mass m q around the chiral limit 
(m q — > 0), which in %PT becomes an expansion in p — 
m„./ A x , the mass of the pseudo-Goldstone boson of sponta- 
neous chiral symmetry breaking over the scale of chiral sym- 
metry breaking A x ~ 4^/^ ps 1 GeV [32]. Because of the 
branch cut in the complex-m 2 plane along the negative real- 
axis, the chiral expansion is not a series expansion (otherwise, 
it would have a zero radius of convergence), but rather an ex- 
pansion in non-integer powers of m% oc m q . 

By writing the dispersion integral as: 



f(ml) = — 




(7) 



it is evident that the second integral can be expanded in integer 
powers of m% /A?, . Hence this term is of analytic form and can 
only affect the values of the LECs. Indeed, the physics above 
the scale A x is not described by %PT and therefore its effect 
should be absorbable by the LECs. 

The second integral generates an infinite number of analytic 
terms, while the number of LECs to a given order of the cal- 
culation is finite. The higher-order analytic terms are present 
and not compensated by the LECs at this order, but their ef- 
fect should not exceed the uncertainty in the calculation due 
to the neglect of all the other higher-order terms. That is, 
the second integral can be dropped, while the resulting cutoff- 
dependence represents the uncertainty due to higher-order ef- 
fects. 

One purpose of imposing a cutoff of order of 1 GeV is 
to investigate the convergence of the expansion without ac- 
tually computing any of the higher-order contributions. This 
is one the main goals of FRR- indeed, motivated by the ab- 
sence of rapid curvature in all hadronic observables at larger 
quark masses, FRR aims to resum the chiral expansion in the 
expectation of improving the convergence of the residual se- 
ries (HI miM HI HI. From the formulae shown in Eqs. <g) 
- © in the previous section, it is clear that the 'sharp cutoff 
FRR is equivalent to the cutoff pion-mass dispersion relation: 
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where n indicates the number of subtractions around the chi- 
ral limit. In this work, the main aim is to see at which values 



IV. NUCLEON PROPERTIES AT O (p 3 ) 

At chiral order p 3 , the imaginary parts of the nucleon mass, 
the proton and neutron AMMs, and the magnetic polarizabil- 
ity of the proton were computed in Ref . [ 1 8 o 
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where Mjv ~ 939 MeV is the physical nucleon mass, 
e 2 /47r ~ 1/137 is the fine-structure constant, and the fol- 
lowing dimensionless variables are introduced: 
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(10) 



The expression for the mass comes from the graph in Fig. |2] 
with leading (pseudo-vector) irNN coupling. The expres- 
sions for the AMMs and polarizability come from graphs 
obtained from Fig. [2] by minimal insertion(s) of 1- and 2- 
photons, respectively. 

The corresponding heavy-baryon expressions at order p 3 
can be obtained by keeping only the leading in 1/M/v term 
(i.e, A ~ y/—r, etc.): 
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1 The original expressions of 11811 contain misprints: Eqs. (10)-(12) miss an 
overall factor of 4, while Eq. (14) misses a factor of r. 



The full, renormalized result for a given quantity is obtained 
by substituting these imaginary parts into the dispersion rela- 
tion of Eq. ([H]). The number of subtractions required in each 

I 



case differ: n — 2 for Mjv, n = 1 for AMMs, and no sub- 
tractions for polarizability. The resulting expressions read as 
follows: 
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The heavy-baryon expressions can be obtained by picking 
out the leading in 1/Af/v term, or equivalently, by substitut- 
ing the corresponding imaginary parts from Eq. (ITTl . into the 
dispersion relation. In the latter case, the same integral is en- 
countered in all of the examples: 



contributions is shown at the physical value of the pion mass: 
m^ ~ 139 MeV. Each quantity (mass, isovector AMM, and 
polarizability) is presented in a separate panel, where the re- 
sults with and without the heavy-baryon expansion are dis- 
played. This figure illustrates the following two features: 
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All of the above quantities to 0(p 3 ) in HB^PT are given by 

this integral, up to an overall constant, and a factor of m 2n . n 

is the number of subtractions (or pertinent LECs) at this order. 

In Fig.[3j the resulting cutoff-dependence of the above loop 
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As seen from Eqs. (fT2l and 1131 . the residual cutoff- 
dependence in HB%PT falls off as 1/A in all of the 
considered examples, while the dependence in the case 
of B^PT behaves as 1/A 2 for M N , and as 1/A 4 for 
both AMMs and (3 p despite the presence of the posi- 
tive powers of A in Eqs. ( fT2b - d). The HB^PT re- 
sults have a stronger cutoff-dependence than the B%PT 
results, indicating a larger impact of the unknown 
high-energy physics to be renormalized by higher-order 
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FIG. 3. (color online). The cutoff-dependence of leading-order loop 
contributions to various nucleon quantities (mass, isovector AMM, 
and proton's magnetic polarizability) calculated in HBxPT (blue 
dashed curves) and BxPT (red solid curves). 



LECs. Hence, the HB results are at bigger risk of pro- 
ducing a large contribution from the high-momentum 
region. Note that, although not necessarily immediately 
apparent from Fig. [5] the cutoff-dependence of the rel- 
ativistically improved chiral formulae can be obtained 
from Eqs. ( fT2h - d), and in the heavy-baryon case, from 
the chiral formulae from the FRR literature; e.g. see 
Refs. J25I2J1. 
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The HB- and BxPT results are guaranteed to be the 
same at small values of A (and not only at A = 0), 
as can be seen by taking derivatives of Eq. ^ with re- 
spect to A 2 , at A = 0. However, at finite values of A 
the differences are appreciable. Observing significant 
differences for A of order m n , as in the case of f3 p , indi- 
cates that the size of the 1/Mjv terms is largely under- 
estimated in HBxPT. 



V. MATCHING: CONFRONTING LATTICE RESULTS 

Eventually, the xPT results must be matched to the underly- 
ing theory - or in practice - fitted to experimental and lattice 
QCD simulation results. In this section, it will be demon- 
strated that chiral extrapolations based on the relativistic ex- 
pressions of Eqs. (fT2h - d) are more stable with respect to 
cutoff variation. 

The case of magnetic polarizability is interesting. Since 
there are no unknown LECs at leading order, it constitutes a 
genuine prediction. Unfortunately, there are no lattice results 
for this quantity, while the experimental data are largely un- 
certain (see Ref. [35] for a recent discussion). One thing that 
the data indisputably show is that /3 P is small compared to the 
electric polarizability, and positive, which seems to be more 
consistent with the HBxPT result. However, it is well known 
that f3 p must have a large positive contribution from the ex- 
citation of the A(1232) resonance [35], which can only be 
accommodated if the chiral loops are negative and partially 
cancel it out. 

For the nucleon mass, one does not expect much difference 
between HB- and BxPT around the physical pion mass, based 
on Fig. [3] For larger pion masses, however, the difference 
becomes significant, and may affect the fit to lattice results as 
is shown in what follows. 

In Figs. |4] and [5] chiral extrapolations of recent lattice re- 
sults from PACS-CS [29] and JLQCD JU are presented. The 
different panels correspond to different values of the cutoff 
A, while the dashed and solid curves correspond to HB- and 
BxPT fit at order p 3 . The values of the fit parameters obtained 
using PACS-CS and JLQCD results are shown in Tables [Jand 
im respectively. 

The PACS-CS results were generated using non- 
perturbatively C(a)-improved Wilson quark action at a 
lattice box size of ~ 2.9 fm, but the set only contains five 
points, and there is a large statistical error in the smallest 
to 2 point. The JLQCD results were generated using overlap 
fermions in Nf = 2 QCD. The lattice box length for each 
simulation result is ~ 1.9 fm, with a corresponding lattice 
spacing is 0.118 fm. The box size is small compared to that 
of the PACS-CS simulations, but the statistical uncertainties 
in each point are also smaller. For simplicity, the fits also 
neglect possible finite-volume corrections. 



Previous studies in HBxPT have indicated that, for lattice results extending outside the chiral power-counting regime, 
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FIG. 4. (color online). Comparison of chiral extrapolations of the nucleon mass for HBxPT (black dashed curves) compared to the relativistic 
formula of BxPT (blue solid curves) from Eq. dl2a) , at various values of sharp-cutoff scale of A. The extrapolation based on PACS-CS results 
(U, box size: 2.9 fm. 
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FIG. 5. (color online). Comparison of chiral extrapolations of the nucleon mass for HBxPT (black dashed curves) compared to the relativistic 
formula of BxPT (blue solid curves) from Eq. dl2at . at various values of sharp-cutoff scale of A. The extrapolation based on JLQCD results 
I3CTI - box size: 1.9 fm. 



the optimal value of the FRR scale A is of the order Awl 
GeV 113311 . Clearly, Fig. [4] shows agreement that the best 
heavy-baryon result is obtained for Awl GeV. The heavy- 
baryon extrapolation is much more sensitive to changes in the 
FRR scale A compared with the BxPT extrapolation, in agree- 
ment with Fig. [3] 

For small values of A, the HBxPT and BxPT results are 
similar, since the chiral loops are suppressed. An almost- 
linear fit eventuates in both cases. This is not ideal, as Fig. [4] 
indicates that neglecting the chiral curvature leads to a poor 
fit of the low-energy lattice results. For larger values of A, the 
heavy-baryon extrapolation struggles to fit the lattice results 
due to large curvature in the heavy pion-mass region. The rel- 
ativistic extrapolation appears to produce a more stable fit to 
the lattice results across a range of values of cutoff scale A. 

The importance of accommodating chiral curvature is even 
greater in the case of observables with lower-order leading 
nonanalytic terms in their chiral expansions, such as the mag- 
netic moment of the nucleon. In the case of the AMMs, recent 
lattice QCD simulations by the QCDSF collaboration are used 
& 



The results from QCDSF were generated using Nf = 2 
and the O (a) -improved Wilson quark action, with box sizes 
ranging from 1.7 to 2.9 fm. To ensure that the lattice results 



TABLE I. A comparison of the values of the fit parameters Mn an d 

o 

Ci for HB^PT and relativistic BxPT for various values of regulariza- 
tion scale A, based on lattice results from PACS-CS. 

A(GeV) M N (GeV) Mjv(GeV) c x (GeV" 1 ) c x (GeV" 1 ) 



0.5 


0.953(22) 


0.963(22) 


-0.709(9) 


-0.660(9) 


1.0 


0.856(22) 


0.896(22) 


-0.970(9) 


-0.840(9) 


3.0 


0.717(22) 


0.833(22) 


-1.278(9) 


-0.981(9) 



from QCDSF give a reasonable approximation to the infinite 
volume limit, the following restrictions are applied: L > 1.5 
fm and m^L > 3. There are nine lattice points that satisfy 
these criteria from the original set of results. Additionally, the 
isovector combination of the nucleon (p — n) is used in lattice 
QCD to avoid calculating the disconnected loops that occur 
in full QCD, which are computationally intensive, since they 
involve the calculation of all-to-all propagators. In general, 
diagrams contributing to the AMM of a hadron include pho- 
tons coupling to sea-quark loops. In the special case of the 
isovector, the diagrams that include these disconnected loops 
cancel. 



♦ Lattice results 
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— FRR rel BxPT: A=0.5 GeV 
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FIG. 6. (color online). Comparison of chiral extrapolations of the anomalous magnetic moment for HBxPT (black dashed curves) compared 
to the relativistic formula of BxPT (blue solid curves), for three different values of A. The extrapolation based on QCDSF results 13 III , box 
size: 1.7 - 2.9 fm. 



TABLE II. A comparison of the values of the fit parameters Mn an d 

o 

C\ for HBxPT and relativistic BxPT for various values of regulariza- 
tion scale A, based on lattice results from JLQCD. 



-ETB" 



~B 



oHB 



A(GeV) M N (GeV) Mjv(GeV) l x (GeV" 1 ) Cj (GeV 1 ) 

0.5 0.986(8) 0.997(8) -0.676(1) -0.627(1) 

1.0 0.880(8) 0.924(8) -0.943(1) -0.811(1) 

3.0 0.721(8) 0.851(8) -1.265(1) -0.959(1) 



TABLE III. A comparison of the values of the fit parameter « isov and 
a,2 for HBxPT and relativistic BxPT for various values of regulariza- 
tion scale A. 



A(GeV) 



7HB - 



a% B (GeV" 



2 ) a? (GeV" 2 ) 



0.5 5.23(5) 4.13(5) 0.70(10) 0.14(10) 

0.8 6.27(5) 4.28(5) 3.42(10) 0.65(10) 

1.0 6.68(5) 4.31(5) 5.46(10) 0.83(10) 



The BxPT integral results from Eqs. (1 12b) & (112c) can be 
adapted for chiral extrapolation of the nucleon isovector by 
taking the difference between the proton and neutron AMM 
formulae. In addition, a linear term in m^ is added with a free 
fit parameter a^, which plays the role of compensating some 
of the high-momentum contributions: 



,(m;,A 2 ) = K p (m;,A^)-K n (m;,A 2 ) + a 2 



(14) 



Without it, the heavy-baryon result is just a straight line, and 
the differences between the two frameworks are irreconcilable 
in this range of pion masses. 

The corresponding chiral behavior of the isovector nucleon 
AMM is shown in Fig. [6] for three different values of A, with 
corresponding values of the fit parameters shown in Table [III] 
Note that the curvature of the extrapolation using a sharp- 
cutoff regulator with A = 1.0 GeV is already large. This 
is a consequence of the leading-order nonanalytic behavior of 
the AMM occurring at a lower chiral order (~ m w ) than for 
the nucleon mass (~ m^). However, the extrapolation using 
the relativistic BxPT formulae of Eqs. ( 112b) & ( 112c) is com- 
paratively insensitive to changes in the FRR scale A. This 
indicates that the BxPT formulae are largely independent of 
the ultraviolet behavior. 

As shown in the previous section the naturalness problem 
doesn't arise in the AMM, in the same drastic way as it does in 
the polarizability. However for a broader range of pion masses 
the problem starts to show. In the case of AMM, a^ B quickly 
becomes larger with increasing A, in order to accommodate 
the large higher-momentum contributions. 



VI. CONCLUSION: WHEN HEAVY-BARYON FAILS 



The HBxPT and BxPT can be viewed as two different ways 
of organizing the chiral EFT expansion in the baryon sector. 
While the heavy-baryon expansion is often considered to be 
more consistent from the power-counting point of view, it ap- 
pears to be less natural. Certain terms that are nominally sup- 
pressed by powers of m^/M^, and hence dropped in HBxPT 
as being 'higher order' , appear to be significant in explicit cal- 
culations. 

The problem is more pronounced in some quantities and 
less in others. To quantify this, one needs to note the power 
of the expansion parameter at which the chiral loops begin 
to contribute to the quantity in question. For the considered 
examples of the nucleon mass, AMMs, and polarizability, this 
power index is 3, 1, and —1, respectively. The smaller the 
index, the greater is the difficulty for HBxPT to describe this 
quantity in a natural way. For quantities with a negative index, 
a dramatic failure of HBxPT is expected. 

The negative index simply means that the chiral expansion 
of that quantity begins at lowest order with negative powers of 
-inn- Apart from polarizabilities, the most notable quantities of 
this kind are the coefficients of the effective-range expansion 
of the nuclear force. As is known, the non-relativistic x?T in 
the two-nucleon sector [36] failed to describe these quantities 
03711 . thus precluding the idea of 'perturbative pions' in this 



sector^ The present work encourages one to think that B^PT 
can solve this problem, as is the case for nucleon polarizabili- 
ties. 

It certainly is important to understand the origin of the ap- 
parent difficulty of the HB expansion in estimating quantum 
corrections in certain observables, at least in a way it has been 
understood for the scalar form factor [7J. Here, it was only 
shown that the behaviour comes from the "soft" momentum 
region, and a criterion for determining the region relevant to 
this problem was conjectured. The understanding of the origin 
presents a challenge for future studies. 



mass dispersion relation (as seen in Eq. (fl3): 



/K) = -i/d^ 



Im/(i) 



TO 



(Al) 



with / being a function of m q . The above allusion to GOR 
implies that this relation is valid for small m q only. However, 
if its validity were assumed for all m q , the issue of its conver- 
gence for a given quantity / could be considered as follows. 

The unsubtracted dispersion relation implies that 
/ ex 1/rriq, for large m q , which obviously cannot be 
true for every quantity. For example, In the case of the 
nucleon mass / = Mm, it is expected that / oc m q for large 
to . Therefore, the relation needs to be subtracted at least 
twice, i.e: 
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M N (m q ) =Mn +aim q / -j 



dt lmM N (t) 



t 2 t 



TO 



(A2) 



with M n and a\ being the subtraction constants. 

In another example, the nucleon AMM should behave a 
constant for large m q and hence one subtraction will suffice: 



K N (m q 



K N 



dt ImKjv(i) 

t t — TO„ 



(A3) 



Appendix A: The quark-mass dispersion relation and 
subtractions 



According to the Gell-Mann-Oakes-Renner relation 



(GOR) [39], to 2 oc TO 9 for a light-quark mass m q . Thus, the 
pion-mass dispersion relation can be translated into a quark- 



where k^ is the AMM in the chiral limit, playing the role of 
the subtractions constant. 

The sufficiency of these subtractions is confirmed in 
leading-order %PT calculations. At higher chiral orders, more 
subtractions are needed as new low-energy constants arise to 
play the role of the subtraction constants. The above analy- 
sis determines only the minimal number of subtractions for a 
given quantity. 
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